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ABSTRACT 

Observational evidence accumulated over the past decade indicates that accretion discs 
in X-ray binaries are viscously stable unless they accrete very close to the Eddington 
limit. This is at odds with the most basic standard accretion disc theory, but could be 
explained by either having the discs to be much cooler whereby they are not radiation 
pressure dominated, or by a more sophisticated viscosity law. Here we argue that 
the latter is taking place in practice, on the basis of a stability analysis that assumes 
that the magneto-rotational-instability (MRI) responsible for generating the turbulent 
stresses inside the discs is also the source for a magnetically dominated corona. We 
show that observations of stable discs in the high/soft states of black hole binaries, 
on the one hand, and of the strongly variable microquasar GRS 1915-1-105 on the 
other, can all be explained if the magnetic turbulent stresses inside the disc scale 
proportionally to the geometric mean of gas and total pressure with a constant of 
proportionality (viscosity parameter) having a value of a few times 10~^. Implications 
for bright AGN are also briefly discussed. 
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1 INTRODUCTION 



According to the widely adopted standard solutions 
(Shakura & Sunyaev 1973), luminous accretion discs close 
to the Eddington rate should be radiation pressure domi- 
nated and therefore unstable to perturbations of both mass 
flow (Lightman & Eardley 1974) and heating rate (Pringle, 
Rees & Pacholczyk 1973) for the commonly adopted assump- 
tion that viscous stresses within the disc are proportional to 
the total (gas-|-radiation) pressure (a- viscosity prescription). 
Taking into account the stabilizing efi'ect of radial advec- 
tion near the Eddington rate, often modeled with the "slim 
disc" solutions (Abramowicz et al. 1988), a limit cycle-type 
of behaviour should be expected, which has been confirmed 
by numerical simulations of time-dependent discs (Taam & 
Lin 1984; Honma, Matsumoto & Kato 1991; Szuszkiewicz & 
Miller 1997). 

These instabilities may also operate in accretion discs of 
supermassive black holes. It is now believed that the main 
phase of growth of these black holes must have occurred 
in short-lived episodes of near-Eddington accretion (Yu & 
Tremaine 2002; Merloni 2004; Hopkins et al. 2006), most 
likely associated with bright quasar phases. The lifetimes 
and duty-cycles of such luminous objects, crucial to under- 
stand their cosmological evolution, may be influenced by the 
radiation-pressure instabilities (or lack thereof) of accretion 
discs. 

Galactic black holes in binary systems provide an im- 
portant laboratory to study these instabilities on humanly 



observable time scales. Among all of them, only one - the 
famous microquasar GRS 1915-1-105 - displays a variabil- 
ity which can be explained by some kind of viscous insta- 
bility (Belloni et al. 1997; Nayakshin, Rappaport & Melia 
2000; Janiuk, Czerny & Siemiginowska 2000; for a recent 
review and a more comprehensive list of references on this 
enigmatic source, see Fender & Belloni 2004). On the other 
hand, the vast majority of transient galactic black holes are 
observed to be stable in a disc-dominated state (the so-called 
high/soft state, or thermal dominant state, see McClintock 
& Remillard 2006) at luminosities which are a few tenths of 
the Eddington one, which is hard to reconcile with the orig- 
inal viscosity prescription of stresses proportional to total 
pressure. 

Two physically plausible ways to make luminous ac- 
cretion disks stable are well known. First of all, the discs 
may be cooling much faster than the standard solution as- 
sumes due to an additional rapid energy transfer from the 
disc mid-plane into a corona, a jet or a wind, so that the 
radiation pressure simply never dominates in the disk (e.g., 
Svensson & Zdziarski 1994). Secondly, the anomalous vis- 
cosity of accretion discs, now understood to be due to MRI, 
may not scale with the total disc pressure (e.g., Lightman & 
Eardley 1974) . Numerical MHD simulations of turbulent ac- 
cretion flows are the most promising tools for differentiating 
between these possibilities from first principles (see, in par- 
ticular, Sano et al. 2004 and references therein). However, 
due to immense numerical challenges, one will have to wait 



2 Merloni & Nayakshin 



until global 3-D radiative simulations are performed over a 
large enough range of parameter space before a clear answer 
will emerge. 

In this paper we investigate the time-dependent evolu- 
tion of magnetized accretion discs by means of semi-analytic 
techniques and discuss their relation to observations. In 
particular, we suggest plausible scalings for the fraction of 
power dissipated in the corona and the viscosity law, and we 
perform a local viscous disc stability analysis, later checked 
with timo-dcpcndcnt disc models. Wc find that the observed 
stability properties of galactic black hole discs are most 
likely explained by a modified viscosity law rather than by 
an extra coronal or jet cooling. 



2 VISCOSITY LAW 

The nature and extent of the posited limit cycle instabilities 
at high accretion rates depend critically on the poorly under- 
stood prescription for the viscous torques: it is well known 
that assuming viscous stresses scale proportionally to the gas 
pressure results in accretion discs which are stable through- 
out, oven at the highest accretion rates (Lightman and Eard- 
ley 1974; Stella and Rosner, 1984). In fact, our ignorance of 
the physical mechanisms giving rise to the disc viscosity, and 
in particular of its exact scaling, has led many authors to 
consider the outcome of radiation pressure dominated discs 
obeying a more general prescription for the viscous stresses 
tr,j, (Taam & Lin 1984; Szuszkiewicz 1990; Honma, Mat- 
sumoto & Kato 1991; Watarai & Mineshige 2003): 



trti 



tot ^gas ) 



(1) 



with ao constant, where 

J'tot = Pgas + Prad = 2p m-'-kT + aT'^l?, (2) 

is the sum of gas plus radiation pressure (for hydrogen rich 

material), T is the mid-plane disc temperature and is 
the proton mass. Within this approach, the parameter /i 
can take any value between (stresses proportional to total 
pressure) and 2 (stresses proportional to gas pressure only) . 



2.1 Magneto-Rotationally Unstable discs and the 
generation of the corona 

Numerical studies of the last decade have shed new light on 
the nature of viscosity in accretion discs, by elucidating the 
crucial role of MHD turbulence for their enhanced transport 
properties. Since magneto-rotational instability (MRI; see 
Balbus & Hawley 1998, and reference therein) is the primary 
driver of the angular momentum transfer in the discs, the 
turbulent magnetic stresses scale with magnetic pressure, 
and therefore tr^, oc Pmag, where Pmag = |Sdisc|^/87r is the 
magnetic pressure inside the disc. 

One of the main open issues in the physics of black 
hole accretion discs is the relationship between the disc 
MRI-driven turbulent viscosity and the generation of the 
hot coronae that are usually postulated in order to explain 
the observed X-ray emission (Liang & Price 1977; Galeev, 
Rosner & Vaiana 1979; Blackman & Field 2000; Kuncic 
& Bicknell 2004). Phenomenological models usually assume 
that at each radius, a fraction / of the internally generated 



power is transferred vertically outside the disc, and pow- 
ers a magnetically dominated corona (Haardt & Maraschi 
1991; Svensson & Zdziaxski 1994). As customary (see e.g. 
Svensson & Zdziarski 1994; Merloni 2003), we assume that 
in MRI-turbulent discs such a fraction / of the binding en- 
ergy is transported from large to small depths by some form 
of collective mean electromagnetic action (Poynting flux). 
One should always keep in mind, however, that this is by no 
means the only way in which energy can be removed non- 
radiatively from the optically thick disc (for an alternative, 
see e.g. Tagger & Pellat 1999). 

We can now estimate the vertical Poynting flux, F-p, in 
the simplest way, assuming that _fp ~ t'D-Pmag, where is 
the upward drift velocity of a magnetic flux tube within the 
disc. In Merloni (2003) it was argued that vd should in gen- 
eral be of the order of the Alfven speed v\. This translates 
into the following expression for the fraction of power dissi- 
pated in the corona, uniquely relating this quantity to the 
magnetic disc viscosity parameter qq (Merloni 2003; Hirose, 
Krolik & Stone 2006): 



(3) 



where /3 = 1/(1-1-^), is the ratio of gas to total pressure, ^ = 
fVad/fgas is the ratio of the radiation to the gas pressure. 
Note that in this approach / is an implicit function of radius, 
through the radial dependence of the pressures. 

Recent progress in numerical studies of the disc-corona 
coupling has been made by simulating a gas-pressure domi- 
nated local patch of an accretion disc (with vertical gravity 
included) in which heating by dissipation of the MHD tur- 
bulence is balanced by radiative cooling (Hirose, Krolik & 
Stone 2006; see also Miller & Stone 2000). In broad accor- 
dance with eq. (3), it was found that the fraction of power 
released outside the disc main body was less than about 
10% for a measured stress parameter of qq « 0.02. How- 
ever, due to the increased magnetic pressure support in the 
upper disc layers, most of the Poynting flux emerging from 
the disc main body is dissipated helow the photosphere, and 
therefore cannot be directly associated with the observed 
hot, optically thin X-ray emitting plasma. Obviously, global 
radiative simulations are needed to assess the role of long- 
wavelength Parker instability modes, and the scaling with 
the radiation pressure predicted by eq. (3) for the generation 
of genuinely hot coronae from disc magentic fields. 



3 STABILITY ANALYSIS 

Let us consider now the general case of the viscosity pre- 
scription (1), with < /i < 2. We can calculate analytically 
the value of the disc parameters at which the instability 
sets in by studying the stability properties of the stationary 
solution. In the one (vertical) zone limit, the equation for 
hydrostatic equilibrium in the vertical direction is 



GM'SH 
2i?a 



(4) 



while the angular momentum conservation equation reads: 



pM/2p(l-M/2) 
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(5) 
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Figure 1. The critical value of the gas to total pressure ratio, 
/3crit I below which the disc is unstable as a function of the viscosity 

parameter ao and index fi. Calculations arc for the case of no large 
scale magnetic field torques, and for a black hole mass A/bh = 
1OM0. The solid lines mark the values /3crit = 0.1,0.2,0.4,0.8. 
The horizontal dashed line mark the // = 1 case discussed in 
section 2. 



where the function J{R) — (1 — ^Rin/R), with Bin — 
3Rs, describes the Newtonian approximation of the no- 
torque at the inner boundary condition for a disc around a 
Schwarzschild black hole. Finally, the energy balance equa- 
tion is given by 

p _cP,ad 3 3nlMJ{R){l - f) 

tt ^ on 

where Frad is the vertical radiation flux and we have / oc 
from eq. (3). 

By differentiating logarithmically the above expres- 
sions, together with the equation of state (2), with respect 
to M, we obtain four equations in terms of the logarith- 
mic derivatives of Ptot, P, T and H, that can be solved for 
dlogp/dlogM and dlog H/dlog M as functions of /3 and /. 
The lower turning point of the Af (S) curve, which indicates 
the local instability condition is then found by setting 

dlogS ^ dlogp dlogH ^ ^ 
dlogM dlogM dlogM 
which in turn gives: 

. _/PgaA _ 7/i(2 - 3/) - 16(1 - /) 
~ I Ro. ) cHt " 7m(2 - 3/) - 40(1 - /) ■ 

The instability sets in at a transition radius which can be 
found by solving the following algebraic equation: 

(o,ro)"»'m'""(l-/)«" (9) 



where we have defined rn = Mbh/Mq. 

In the limit of / = 0, we obtain the known result (see 
Szuszkiewicz 1990) that, for the instability condition to be 
satisfied, the gas pressure needs to be smaller than 0.4 times 
the total pressure if /i = 0, and just 1/13 times if /i = 
1. Thus, the larger is yu, the more difficult is to excite the 
instability. By using eq. (3) to relate / and /?, we can then 
explore the whole space spanned by the parameters ao and 
fj.. In Figure 1, where we plot as contours the surface of 
critical gas to total pressure ratio Pait, while in Figure 2 
we plot the corresponding critical values of the accretion 
rate, defined as m s eM(? / L^m, where e is the radiative 
efficiency"'^. 

First of all, as it is well known, there exists a region of 
parameter space for which the disc is always stable. Such 
a region correspond to the upper left corners of figures 1 
and 2, i.e. for large values of /u, implying stresses coupled 
more with the gas, rather than the total pressure, and small 
magnetic viscosity parameters qq. 

On he other hand, if viscous stresses are proportional 
to total pressure (/it = 0), the presence of a corona as a 
sink of energy has a stabilizing effect on the disc (Svensson 
& Zdziarski 1994), and the critical accretion rate scales as 
rhcrit oc (1 — /)~^^*. However, if ao is large enough the disc 
can be unstable even for large values of ji (upper right corner 
of figures 1 and 2). 

3.1 Observational constraints 

It is an observational fact that only the microquasar GRS 
1915-1-105, the most luminous transient black hole binary in 
our galaxy, displays strong variability on timescales compat- 
ible with the limit cycle behaviour expected from thermal 
instabilities (see Belloni et al. 1997; Nayakshin et al. 2000; 
Janiuk et al. 2000,2002). Many other systems are known to 
be stable up to luminosities which are a few tens of per- 
cent the Eddington luminosity. This is clearly impossible 
for standard viscosity prescriptions (/i = 0), unless the vis- 
cosity parameter ao (and thus the fraction of power released 
into the corona) is high enough. But, as figures 1 and 2 show, 
there is only a very limited region of the parameter space for 
which this is possible, and within such a region, the fraction 
of power released outside the optically thick disc is always 
large, in contradiction with the observed spectral properties 
of high/soft state black hole binaries. 

On the other hand, if /i ?a 1, as argued, for example 
in Merloni & Fabian (2002) or Merloni (2003), we can have 
accretion discs which are stable at more than half of the Ed- 
dington luminosity, and this holds over a very large range in 
possible viscosity parameters ao - This can indeed explain the 
observed lack of unstable, disc-dominated galactic black hole 
binaries up to luminosities of about a half of the Eddington 
one (Gierlihski & Done 2004), while at the same time allow- 
ing for some instability, taking in place only in those sys- 
tems which are constantly accreting very close to (or above) 

^ In the following, for the sake of consistency with our Newtonian 
inner boundary condition, we will assume e = eo = 1/12. How- 
ever, it should be kept in mind that for a black hole of arbitrary 
spin, the true critical value of L/L-g,dd should be rescaled by a 
factor e/eo with respect to what shown in Fig. 2. 
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Figure 2. The critical accretion rate in units of Eddington, 
riicriti calculated for a radiative efRciency of eg = 1/12, above 
which a disc is unstable is plotted as contours as a function of 
the viscosity parameter qo and index fi. Calculations arc for a 
black hole mass A/bh = IOMq. The solid lines mark the values 
"icrit = 0.15,0.2,0.3,0.5,0.9. The horizontal dashed line mark 
the jU = 1 case discussed in section 2. 



the Eddington rate, as supposedly is GRS 1915+105 (Done, 
Wardzinski & Gierliriski 2004). 

Wc can extend this result to the case of supermassive 
black holes (SMBH). It is well known that for any value 
of /X and ao the critical accretion rate scales with mass as 
7V^-i/8^ Then for a 10^ solar masses black hole the 
limit cycle instability should set in at a luminosity, in units of 
Eddington, which is just one tenth of that of a binary black 
hole, both viscosity law and black hole spin being the same. 
If indeed most of the bright, high redshift Quasars that dom- 
inates the growth history of SMBH have L/L^dd > 0.1, 
as recently proposed (Vestergaard 2004; McLure & Dunlop 
2004), then they should be undergoing limit cycle instability. 



4 TIME DEPENDENT EVOLUTION 

In order to test the above conclusions based on the stabil- 
ity analysis of stationary solutions, we have carried out a 
set of time dependent simulations for different values of the 
parameters qq and rh. From the equation of conservation of 
mass and angular momentum, we can write (Pringle, 1981; 
Livio & Pringle 1992) 



as _ 3__9_ 

dt ^ RdR 



R 



(10) 



The time-dependent energy equation needed to close 
the set must take into account large radial gradients of 
temperature, and therefore includes a number of additional 
terms with respect to the standard heating and cooling 
terms of the stationary solution. The form of the energy 



equation that we use follows the formalism of Nayakshin et 
al. (2000): 
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where 7 is the ratio of specific heats (7 = 5/3) and Fs is 
given in Abramowicz ct al. (1995). Here, the radial velocity 
vr induced by viscous stresses is given by (Eq. 5.7 of Frank 
et al. 1992): 



Vr = 



_d_ 

OR 



1/2] 



(12) 



The terms on the left hand side of Equation (11) repre- 
sent the full time derivative (e.g., d/dt + VRd/dR) of the gas 
entropy, while the terms on the right are the viscous heat- 
ing, the energy flux in the vertical direction, the diffusion of 
radiation in the radial direction, and the viscous diffusion 
of thermal energy. Following Cannizzo (1993; and references 
therein), we take J = 2cpu{T,/R)[d{RdT/dR)/dR] to be 
the radial energy flux carried by viscous thermal diffusion, 
where Cp is the specific heat at constant pressure. F+ is the 
accretion disk heating rate per unit area, and is given by 



F+ = (9/4)!.Sn|,(l-/). 

The radiation flux in the radial direction is 



Tt 



dR 



(13) 



(14) 



where tt is the optical depth of the disk, tt = k;E/2, and k 
is the radiative opacity (assumed here to be dominated by 
electron scattering opacity). Finally, the radiative cooling 
rate in the vertical direction is given by 



(15) 



The results of our time-dependent simulations for the 
case n = 1 and Mbh = 14Mo (as appropriate for GRS 
1915-1-105) are shown in figure 3 as a set of lightcurves plot- 
ted over a time typically of the order of the viscous time at 
the outer domain boundary. Shown separately are the total 
disc (thin red solid lines) and coronal (blue dashed lines) 
emissions, together with their sum (thick black solid lines), 
representing the total dissipated energy in the disc-corona 
system. These results broadly confirm the analysis presented 
in section 3, in that we observe stable discs at m > 0.4 and 
0.2 for ao = 0.03 and 0.1, respectively (see figure 2). 

We find in general a relationship between the ampli- 
tude of the instability and its duty-cycle (i.e. the ratio of the 
"burst duration" to the period of the oscillation), whereas 
small duty cycles are associated with large amplitude vari- 
ability. The amplitude itself is much smaller than in stan- 
dard a-discs, and grows with the ratio rh/rhcrit- For accre- 
tion rates just above the critical values, the large cycles im- 
ply a "refilling" time of the inner disc faster than the viscous 
time at the transition radius. In fact, we found that the vis- 
cosity prescription, /x w 1, reproduces many features of the 
phenomenological one introduced by Nayakshin et al. (2000) 
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Figure 3. Set of lightcurves for magnetized discs with tiie vis- 
cosity law corresponding to ^ = 1, calculated for different ac- 
cretion rates (m = 0.2, 0.4, 0.5) and viscosity parameters (oq = 
0.03,0.1,0.2). Each panel shows one lightcurve over a timescale 
approximately equal to one tenth of the viscous time at the outer 
boundary of our discs (r = 200). Thick black lines represent the 
total (disc plus corona) emission, thin red lines the optically thick 
disc emission, and thin blue dashed ones the coronal (X-ray) emis- 
sion. The black hole mass is fixed to IOMq . 



in order to explain the time variability of GRS 1915-1-105 (a 
task which is beyond the aim of this work). 

Another general property of the simulated lightcurves 
is that, thanks to the nature of the magnetic viscosity law 
adopted, the fraction / of power dissipated outside the op- 
tically thick disc is higher when gas pressure dominates, i.e. 
when the inner disc is denser and cooler, thus making coro- 
nal X-ray emission more prominent during the low luminos- 
ity parts of the instability cycles. 

Figure 4 shows a set of simulated lightcurves with qq = 
0.03 and /i = 1 for supermassive black holes of different 
masses. As we noted in section 3, for larger black hole masses 
the instability sets in at lower L/L-Edd- As the duty cycle 
also gets smaller as the instability amplitude grows, then we 
must conclude that the higher the external accretion rate, 
M oc Mm, the more time will be spent by a source in the 
dense, cold, low-luminosity portion of the instability cycles, 
close to its critical luminosity. 



Figure 4. Set of lightcurves calculated with = 1, ao = 0.03, for 
different accretion rates (m = 0.05, 0.1, 0.2) and black hole masses 
(Mbh = 10*^, lO'', 10* Mq). Each panel shows one lightcurve over 
a timescale approximately equal to half of the viscous time at the 
outer boundary of our discs (r = 200). Linestyles are as in Fig. 3. 



5 DISCUSSION 

The stability analysis we have presented here can be used 
to explain the observed properties of accreting stellar mass 
black holes in different states. First of all, a major constraint 
can be derived from the absence of any limit-cycle instability 
in disc-dominated spectra. Those can be effectively selected 
in two complementary ways: by detail modeling, choosing 
those spectra in which the power-law component contribute 
to less than a fixed amount (say, 10%); and by studying 
the luminosity-temperature relation of the disc component, 



selecting those spectra for which Ldis 



T as expected 



from an optically thick disc extending down to the inner- 
most stable orbit (Kubota & Makishima 2004). Modulo the 
distance uncertainty, these studies unambiguously demon- 
strated that accretion discs can be viscously stable up to at 
least idisc/^Edd > 0.3 (Gierliriski & Done 2004). 

From Figure 2, then, we conclude that only two disjoint 
regions of parameter space are consistent with the observed 
phenomenology of black hole binaries. The first corresponds 
to Qfo ^ 0.3 and, within our theoretical framework, it corre- 
sponds to accretion discs sandwiched by powerful, stabilizing 
coronae as postulated by Svensson & Zdziarski (1994). The 
second acceptable region of the parameter space is roughly 
delimited by /i>0.8 and ao <0.15 and would indicate that 
the stability properties of high/soft state black holes are dic- 
tated by a modified viscosity law. 
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The former constraint, however, implying a large frac- 
tion of power transported vertically by magnetic fields, 
seems inconsistent with the spectral properties of the ob- 
served black holes in the high/soft state (but see however 
Hirosc et al. 2006 or Blacs ot al. 2006). Moreover, for the only 
source that exhibits limit-cycle instabilities, the microquasar 
GRS 1915-1-105, one can associate the observed variability 
time-scales with the disc viscous time (Belloni et al. 1997; 
Nayakshin et al. 2000; Janiuk et al. 2000). These models con- 
sistently suggest that ao ~ a few percent. If this is indeed 
a general properties of all accretion discs, a coherent pic- 
ture emerges for black hole binaries in the high/soft state in 
which MRI-drivcn turbulence generates stresses which scale 
with the disc pressure as the geometric mean of gas and 
total pressure (Merloni 2003), i.e. fi ~ 1, with a constant 
of proportionality ao of about a few percent. This in turn 
implies quite low values of / and consequent small contribu- 
tion from the power-law component to the observed X-ray 
spectra. 

5.1 The low/hard state 

The scenario we have outlined above applies only to the 
high/soft state of black hole binaries^. It seems unlikely that 
the soft-to-hard transition can be solely due to a decrease 
of the external accretion rate and a consequent increase of 
the gas pressure-dominated part of the disc, accompanied 
by an increase of the strength of the corona (Merloni & 
Fabian 2002). In fact, low values of ao in the high/soft state 
also imply low values of / ^ V2ao in the low state, if the 
disc physics is not dramatically altered by the m variation. 
Thus, there must be an additional physical mechanism, not 
included in the simple treatment of the coupled magnetized 
disc-|-corona system presented here, which is responsible for 
the observed state change. Obvious candidates are: a ra- 
dial transition to an inner optically thin, radiativcly ineffi- 
cient flow, due to some kind of disc evaporation (Meyer & 
Meyer-Hofmeister 1994; Dullemond & Turolla 1998; Spruit 
& Deufel 2002), or a global rearrangement of a large scale 
poloidal magnetic field, so that the accretion energy can be 
dissipated almost entirely into the bulk flow of a relativistic 
jet (Livio, Pringle & King 2003). 

5.2 Variability in AGN 

Since the radiation-pressure dominated regime sets in at 
smaller dimensionless accretion rates for higher black hole 
masses, large scale variability due to radiation-pressure 
driven viscous instabilities must be more widespread in AGN 
than in stellar mass accreting black holes. However, observ- 
ing such variability in practice is only feasible for lower mass 
super-massive black holes. Indeed, the viscous timescalc, 
t^i^^ = {l/Q){l/aP''^'^){R/H)'\ calculated at the transition 
radius (see eq. 9) is approximately given by 

<visc(nO « 40a-'/'m'/'m'/'f3l°/,'-^^'{l - /3„it)i°/^ s.(16) 

This ranges from a few tens to a few hundreds of years if 
m=10'^, consistent with estimates of switch on/off times for 

^ For a general discussion of black hole binaries low/hard state, 
the reader is referred to MarkoflF (2005) and Narayan (2005). 



some "changing look" AGN (Guainazzi et al. 2005), to 10''- 
10^ years for a 10^ solar masses black hole. This, interest- 
ingly, coincides with the proposed intermittency time of ra^- 
dio galaxies, based on radio sources number counts and on 
the properties of radio galaxies size distribution (Reynolds 
& Begelman 1997). 

Another interesting aspect of the predicted variability 
for AGN concerns their feedback on galaxy formation. In- 
deed, in our model, black holes accreting at a time-averaged 
sub-Eddington rate may spend a fraction of their time ac- 
creting at above-Eddington accretion rate. If AGN feedback 
(for example in the form of powerful relativistic jets) is sig- 
nificant when L Si iEdd, as suggested by the phenomenology 
of microquasars (see Fender, Belloni & Gallo 2004), then 
such a source would produce feedback while on the "hot" 
branch of the S-curve, whereas a completely stable source 
at the same time-averaged accretion rate would not. 



6 CONCLUSIONS 

In this paper we discussed the implications of the observed 
long term stability of the vast majority of black hole X-ray 
binaries in the so-called high/soft state (or thermal domi- 
nant state). This fact, together with the properties of the 
lightcurves of GRS 1915+105, the most luminous galactic 
black hole and the oidy such system to display limit-cycle- 
type of instabilities, can be used to put constraints on the 
nature of the viscosity law in accretion discs. 

We performed a stability analysis for accretion flows 
with a flexible viscosity prescription under the hypothesis 
that the Magneto- Rotational-Instability responsible for gen- 
erating the turbulent stresses inside the discs is also the 
source for a magnetically dominated corona. By varying the 
scaling index of the viscosity law, /j., and its overall normal- 
ization, Qo, we have identified those regions of the parameter 
space for which limit-cycle instability can develop as a func- 
tion of the accretion rate. From this analysis a coherent pic- 
ture emerges for black hole binaries in the high/soft state in 
which MRI-driven turbulence generates stresses which scale 
with the disc pressure as the geometric mean of gas and total 
pressure, i.e. A* ~ 1, with a constant of proportionality ao of 
about a few percent. This in turn implies quite low values 
of / and consequent small contribution from the power-law 
component to the observed X-ray spectra. Provided that the 
viscosity law itself does not change dramatically at lower 
accretion rates, a consequence of this result is that the tran- 
sition to the low/hard state must be caused by additional 
physics, such as the evaporation of, or a global large scale 
magnetic field re- arrangement in the inner portions of the 
disc. 

The scaling with black hole mass of the accretion equa- 
tions is such that limit cycle instabilities should play a role in 
all bright AGN and Quasars. The typical timescalc for these 
oscillations, however, grows with black hole mass faster than 
linearly. Only smaller mass black holes (less than a few times 
10^ solar masses) may have limit cycle instability timescales 
of just a few years. Even considering the uncertainties in the 
absolute values of the predicted critical luminosity that come 
from uncertainties in the black hole mass and spin, it would 
be interesting to search for any evidence of the predicted 
behaviour by looking at distribution of observed disk lumi- 
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nosities in large samples of AGN with available estimates of 
the central black hole mass. 
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